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ABSTRACT 

The Suzaku X-ray satellite observed the young stellar object VI 647 Ori on 
2008 October 8 during the new mass accretion outburst reported in August 2008. 
During the 87 ksec observation with a net exposure of 40 ks. V1647 Ori showed 
a high level of X-ray emission with a gradual decrease in flux by a factor of 
5 and then displayed an abrupt flux increase by an order of magnitude. Such 
enhanced X-ray variability was also seen in X ’’.MM- Newton observations in 2004 
and 2005 during the 2003—2005 outburst, but has rarely been observed for other 
young stellar objects. The spectrum clearly displays emission from Helium-like 
iron, which is a signature of hot plasma (kT ~5 keV). It also shows a fluorescent 
iron Ka line with a remarkably large equivalent width of ~600 eV. Such a large 
equivalent width indicates that a part of the incident X-ray emission that irradi- 
ates the circumstellar material and/or the stellar surface is hidden from our line 
of sight. XMM-Newton spectra during the 2003-2005 outburst did not show a 
strong fluorescent iron Kct line, so that the structure of the circumstellar gas very 
close to the stellar core that absorbs and re-emits X-ray emission from the central 
object may have changed in between 2005 and 2008. This phenomenon may Ire 
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related to changes in the infrared morphology of McNeil’s nebula between 2004 
and 2008. 


Subject headings: stars: formation — stars: individual (V1647 Ori) — stars: 
pre-main sequence — X-rays: stars 


1. Introduction 

Certain young stars dramatically increase their mass accretion rates by orders of mag- 
nitude (from 10~ 7 M o yr -1 to 1O~ 4 M 0 yr -1 ). These events are possibly triggered by thermal 
disk instabilities and are traced by dramatic increases in optical/IR luminosities; Such erup- 
tive pre-main sequence (PMS) stars are crudely classified as either FU Ori (FUor) or EX 
Lupi (EXor) types; the former are characterized by outbursts lasting decades, whereas the 
latter generally display shorter outbursts, of duration a few months or years (see reviews in 
Hartmann & Kenyon 1996). Although PMS stars are notoriously luminous X-ray sources, 
FUor or EXor events have been reported only rarely, such that the level of PMS X-ray 
activity during an outburst is poorly established. 

The young stellar object (YSO) V1647 Ori, which is deeply embedded in the L1630 
dark cloud (d ~400 pc, Anthony-Twarog 1982), underwent a strong opticai/NIR outburst in 
2003 December. This eruption afforded the first opportunity to record the sustained X-ray 
outburst of a rapidly accreting YSO (Kastner et al. 2004; Grosso et al. 2005; Grosso 2006; 
Kastner et al. 2006). Multiple Chandra and XMM-Newton observations through this outburst 
demonstrated that the average X-ra.y flux varied in the same way and on approximately the 
same timescale as the optical and IR brightnesses, although the X-ray flux also varied strongly 
(by a factor of up to 20) on timescales of less than a day. During outburst, the V1647 Ori X- 
ray spectrum had the characteristics of deeply embedded hot plasma (Ah ~4.1 xlO 22 cm -2 , 
kT ~4.2 keV). This result appears to be explained best as star-disk magnetic reconnection 
activity generated in association with the episode of very rapid mass infall. 

Following the VI 647 Ori event, a, couple of YSO systems have experienced mass accretion 
outbursts (V1118 Ori: Audard et al. 2005. Z CMa: Stelzer et al. 2009). Follow-up X-ray 
observations, however, detected little enhancement in the X-ray flux. Skinner et al. (2006. 
2009) reported that two FUor systems with ongoing historical outbursts, FU Ori and V1735 
Cyg. showed very hard spectra of kT > 5 keV plasmas or equivalent despite relatively steady 
X-ray fluxes. Systematic surveys of low-mass PMS (T Tauri) stars in the Orion nebula 
and the Taurus dark cloud (e.g. Preibisch et al. 2005; Telleschi et al. 2007) indicate that 
accreting T Tauri stars are less X-ray active (by a factor of ~2— 3 on average) than non- 
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accreting T Tauri stars, suggesting that mass accretion activity may in fact suppress X-ray 
activity. Hence, the relationship between PMS mass accretion and high energy radiation 
remains unclear. 

V1647 Ori began a new optical/NIR outburst at a certain time between Jan. 2 and Aug. 
26, i.e. , when it was not observable from ground based telescopes due to its conjunction 
with the Sun (Itagaki et al. 2008; Aspin 2008; Aspin et, al. 2009). Our team triggered 
an anticipated Target of Opportunity (ToO) observation of V1647 Ori with Chandra in 
September 18 for 20 ksec, in which the X-ray count rate was double the level observed by 
Chandra on 2004 March 7 (Kastner et al. 2004) during the previous outburst in 2003—2005 
(Weintraub et al. in preparation). This new X-ray eruption offers us the first (and possibly 
unique) opportunity to measure X-ray emission during two outbursts from the same YSO. 
We therefore proposed a ToO observation with the Suzaku X-ray observatory (Mitsuda et al. 
2007), so as to obtain detailed X-ray spectral diagnostics and a ~da,y timescale light curve 
of VI 647 Ori during its latest accretion outburst. With this information, especially derived 
from hard X-ray data above ~2 keV, we can derive the hottest plasma temperature and the 
intrinsic activity variation of V1647 Ori, unaffected by circumstellar absorption; we can also 
constrain models of the plasma geometry and of the response of neutral disk material to 
X-ray irradiation. 


2. Observation 

Suzaku observed V1647 Ori for ~87 ksec on 2008 October 8, when the YSO’s opti- 
cal/NIR fluxes leveled off at a brightness similar to that observed in 2004 Feb. -Mar. (Ojha 
et al. 2008; Aspin & Reipurth 2009). At the time of this observation, the observatory op- 
erated with three sets of telescope and detector systems — the thin-foil X-Ray Telescope 
(XRT) module (Serlemitsos et al. 2007) with the X-ray CCD detector (XIS; X-ray Imaging 
Spectrometer, Koyama et al. 2007) on its focal plane, and the Hard X-ray Detector (HXD: 
Takahashi et al. 2007; Kokubun et al. 2007) with narrow collimator arrays. The XRT+XIS 
system focuses X-rays with a half-power diameter (HPD) of ~2' and covers the X-ray energy 
range up to ~10 keV. The XIS1 camera employs back-illuminated (BI) CCDs for the soft 
band sensitivity down to ~0.3 keV, while the other XIS cameras (XIS0, XIS3) use front-side 
illuminated (FI) CCDs with sensitivity down to ~0.4 keV. The HXD detects hard X-rays 
above 15 keV. 

We started analysis from the cleaned event data in the distribution package, pre- 
processed with the version 2.2.11.22. The cleaned . data excluded events recorded during 
high background periods, mainly when the satellite passed the South Atlantic Anomaly 
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(SAA) 9 times during this observation, and earth occultations, each of which happened for 
~2.1 ksec in every Suzaku orbit (~90 min). The data did not suffer telemetry saturation. 
We thus selected events in good time intervals of all the XISO, 1 and 3 sensors. After this 
selection, the net exposure was 40,446 sec, i.e., ~46% of the observing time. For the XIS 
arf generation, we used the CALDB version released on 2009-02-03, to better reproduce the 
recent growth of contamination on the XISs, in particular XISO. 

We generated time-averaged HXD spectra, following the standard analysis procedure 1 
and using the ver 2.0 tuned background 2 . The spectrum was consistent with the typical 
cosmic X-ray background spectrum and showed no hint of signal from V1647 Ori. 


3. Image Analysis 

Figure 1 shows an XIS color image of the V1647 field, combining all the XIS detectors 
(0+1+3). The XISs detected three bright peaks at the center, top, and bottom of the image 
and four faint sources at ~1.5' north-east, 4' south-east, 4' south, and 8' south-west from 
the central source. The bright sources to the north and south and the faint source to the 
south-east each have only one corresponding X-ray source in the 2002 Chandra observation, 
respectively (source 14, source 8, and source 23 in Simon et al. 2004) within the Suzaku 
pointing uncertainty (<20") 3 . We measured their source positions from a 0.5-10 keV XIS 
image smoothed with the Gaussian function with cr= 2 image pixel (16.7") by weighting- 
photon counts within a 3x3 pixel region centered at the peak pixel. These positions were 
offset, on average, by (AR.A., ADec.) — (—8.3", —2.2") from the coordinates in Simon et al. 
(2004), with a standard deviation of 2.0". We shifted the Suzaku image frame to match the 
Chandra coordinates. 

After these corrections, the absolute coordinates of the central bright source ( 0 : 2000 . £ 2000 ) 
= (5 /l 46 m 13.35 s , -00 w 06'5.7") correspond to those of V1647 Ori within the 3cr error range. 
Faint Chandra sources around V1647 Ori (9, 13, 15) are on the outskirts of the PSF and do 
not show any apparent local peaks at their locations in the Suzaku. image. Their contamina- 
tion of VI 647 Ori in the Suzaku data would therefore appear to be minimal. The soft source 
to the north-east from V1647 Ori is probably a blend of sources 20 and 21 in Simon et al. 
(2004). while the very hard source to the south is a blend of sources 7. 10, and 11, which are 


1 http://heasarc.gsfc. 11 asa.gov/docs/suzaku/analysis/abcl 

“ftp://legacy.gsfc.iiasa.gov/suzaku/clata/background/pinnxb-ver2.0-tuned/ 

Isuzakumemo- 2007-04 — http://www.astro.isas.ac.jp/suzaku/doc/ 
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members of SSV 63, a cluster of class I protostars. The source to the west, which is outside 
of the Chandra FOV. may be associated with 2 MASS J05454282-000T261. 


4. Timing and Spectral Analysis 

The north-east source heavily contaminated the V1647 Ori region below ~1 keV, but its 
contamination was negligible above 4 keV. To minimize the contamination and maximize the 
photon statistics, we defined a source region with a 1.5' radius circle, excluding the circular 
region within 40" of the north-east source (see the right panel of Figure 1). To cancel out the 
remaining contamination, we defined a background region that is symmetric to the source 
region with respect to the north-east source. The background-subtracted light curves and 
spectra presented in Figure 2—4 were generated from the events within these regions. 

Figure 2 shows the background subtracted X-ray light curve of V1647 Ori covering the 
energy range 1-8 keV, combining all the XIS (0+1+3) data, and using a bin size of 2 ks. 
The few' missing data points are time periods where few good-time intervals (<1%) were 
available. The count rate varied strongly, by a factor of >18, during the observation (count 
rate of the maximum bin: ~0.11 counts s _1 , lcr upper limit count rate of the minimum bin: 
~6xl0 -3 counts s -1 ). The X-ray light curve shows a gradual decrease, to a level consistent 
with zero flux, over the first ~60 ks of the observation, then a sharp increase in flux at 
~60 ksec followed by a decay with significant spikes and dips. 

We produced background subtracted soft band (1-4 keV) and hard band (4-8 keV) 
light curves and calculated the hardness ratio (HR) defined by (H — S)/(H + S), where S 
and H are the soft and hard band count rates, respectively (bottom panel of Figure 2). The 
HR was ~ —0.5 between 10 ksec and 60 ksec into the observation and increased to ~0.0 
after ~60 ksec, remaining at that level for the rest of the observation duration. The abrupt 
flux increase at ~60 ksec, accompanying an HR increase, may be explained by an increase in 
kT (see also spectral analysis in the latter paragraphs), while fluctuations in flux after that, 
which are apparently uncorrelated with the HR, caused by changes in the plasma emission 
measure. Such flux and HR variations were seen in the XMM-Newton observation in 2004, 
as well (Grosso et al. 2005). 

We merged spectra of the FI sensors and produced time-averaged FI and BI spectra 
(Figure 3). All the spectra in our analysis were binned so as to obtain at least 20 photon 
counts per bin. We ignored energy bins below' 1.4 keV, where contamination from the north- 
east source exceeds the net signal from VI 647 Ori. The resulting spectra reveal a strong 
fluorescent iron line at 6.4 keV from quasi-neutral iron atoms, as well as the helium-like iron 
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line at 6.7 keV, which is a signature of hot plasma. These features can be seen clearly in 
the Fe Ka line region of the merged XIS spectra (right panel of Figure 3). The spike at 
5.9 keV possibly originates from contamination of the 55 Fe calibration source scattering over 
the CCD chips. 

The spectrum is well fit by a model consisting of a single-component, optically thin- 
thermal plasma emission (APEC) combined with a narrow Gaussian line with the center 
energy fixed at 6.4 keV for the iron fluorescent emission, suffering photoelectric absorption 
by neutral gas (wabs). The best-fit, parameters are listed in Table 1 (“whole”). The best-fit 
plasma temperature (~4.1 keV), elemental abundance (~0.50 solar), hydrogen column den- 
sity (~3.6xl0 22 cm -2 ), and the average observed flux (~4.4 xlO -13 ergs cm -2 s -1 between 
0.5-8 keV) are similar to those inferred for V1647 Ori during the XMM-Newton observation 
in 2004 (Grosso et al. 2005). In the 2.8— 8.0 keV band, the absorption corrected X-ray lu- 
minosity during the Suzaku observation (log L x ~30.9 ergs s -1 ) was twice that seen during 
the XMM-Newton observation (log L x ~30.5 ergs s -1 ). One remarkable difference, however, 
is the firm detection of the iron fluorescent line at 6.4 keV. The equivalent width (EW) in 
2008, ~600 eV, was about a factor of 6 higher than that found from the marginal detection 
of the iron fluorescent line during the XMM-Newton observation in 2004. 

Because the spectrum changed so dramatically near the midpoint of the Suzaku expo- 
sure, we divided the observation into two phases, before and after 56 ksec from the obser- 
vation start, as shown at the top of Figure 2. We fitted these spectra independently with 
the same spectral model used to fit the spectrum extracted for the full 87 ksec. The best-fit 
results are shown in Table 1 (“early” and “late”) and Figure 4. The plasma temperature was 
significantly lower in the early phase, while the other parameters, especially the hydrogen 
column density, do not appear to vary significantly. This would suggest that the plasma 
heated up and/or that a very hot plasma emerged at ~56 ksec after the start of the obser- 
vation. The best-fit EW of the iron fluorescent line was huge, close to 1.5 keV, during the 
“early” phase though the error was large enough to include the best-fit EW in the “late” 
phase. Such large fluorescent iron line EWs have been observed from spatially extended X- 
ray reflection sources that do not contaminate emission from their irradiating sources, such 
as molecular clouds that reradiates X-rays from the Galactic center or a bipolar nebula that 
reradiates X-rays from a massive stellar binary system (e.g.. Koyama et al. 1996: Corcoran 
et al. 2004). 



5. Discussion 


The Suzaku observation was performed between ~40 and ~280 days after the onset of 
the new accretion outburst from V1647 Ori. The first (2004) XMM-Newton observation of 
V1647 Ori during its 2003-2005 outburst was performed with similar timing, ~4 months 
after outburst onset, with net exposure time (37—39 ksec) similar to that of the Suzaku ob- 
servation. Since these two observatories have similar effective areas in the 0.5—10 keV energy 
range, these observations are well suited to comparing the X-ray properties of V1647 Ori 
between the 2003 and 2008 outbursts. 

The flux of X-ray emission from V1647 Ori strongly varied, by at least a factor of 18, dur- 
ing the Suzaku observation. The flux bottom is not well constrained given the contamination 
from the north-east soft source, but the XMM-Newton observation in 2004 showed a similar 
range of variation of between 0.5— 12x 10 -13 ergs cm -2 s -1 . Both observations showed similar 
abrupt flux increases with correspondingly HR increases and no significant HR. variations 
after the increases. The plasma parameters derived from the Suzaku spectra were similar to 
the best-fit result of the XMM-Newton spectrum ( kT ~3.0 keV, jV H ~2.9xl0 22 cm -2 and 
Z ~0.8 solar, model #1 in Table 1 of Grosso et al. 2005). These similarities suggest that 
the new outburst in 2008 was driven by a mechanism similar to that of the outburst which 
started in 2003 October and lasted for ~2 years. Similarities between the 2003—2005 and 
latest outbursts are also seen at near-infrared wavelengths (Aspin et al. 2009). 

While V1647 Ori displayed strikingly similar X-ray activity patterns during its two mass 
accretion outbursts, there appears to be a very wide range of X-ray behavior and pre/post- 
outburst X-ray luminosities among the other FUor and EXor systems observed in X-rays 
at the onset of and/or during eruptions. The FUor systems FU Ori and V1735 Ori emit 
strong X-rays at levels of log L x ~31 ergs s _1 from hot plasma, kT > 5 keV (Skinner et al. 
2006, 2009), similar to that seen from V1647 Ori. On the other hand, several FUor and 
EXor systems did not show significant X-ray activity associated with their mass accretion 
outbursts (EXor — V1118 Ori: Audard et al. 2005, FUor — V1057 Cyg and V1515 Cyg: 
Skinner et al. 2009, Z CMa: St.elzer et al. 2009). These particular X-ray characteristics 
are apparently not directly related to the optical outburst criteria that are used to classify 
objects as FUors or EXors (e.g., outburst duration). The X-ray behavior during outburst 
therefore might be more reasonably connected to the other YSO characteristics, such as 
rotational velocity and/or strength of the magnetic field. 

The Suzaku spectra in 2008 showed an impressively strong fluorescent iron Ka line 
feature at 6.4 keV. Such a large E\Y value has been reported only from a couple of young 
stellar objects, including an extremely young protostar ( X E or IRS 7B, Hamaguchi et al. 
2005) and a YSO undergoing a flare (V1486 Ori, Czesla k Schmitt 2007). The observed 



iron fluorescent EWs (Table 1: “whole” and “late”), are ~10 times as large as that for a 4 
keV plasma irradiating an infinite plane of solar iron abundance in a Monte-Carlo simulation 
(~60 eV: flare height of 0.00 R* in Drake et al. 2008). Changing the plasma geometry - 
increasing the plasma height h and/or viewing angle 0 — only decreases the EW ratio. 
Therefore, special physical conditions are required to explain such a large EW. Drake et al. 
(2008) listed four plausible conditions that can elevate the EWs — (1): supersolar iron 
abundance, (2): disk flaring, (3): line-of-sight obscuration of the irradiating plasma, (4): 
excitation by non-thermal electrons. Among them, the condition (4) would be negligible, 
while the condition (2) can only increase the EW by a factor of 2. The condition (1) is 
unlikely because, in this condition, the photosphere and/or disk abundances would have had 
to increase by a factor of six between 2004 and 2008. Thus, the condition (3) appears to be 
the most likely explanation. 

Assuming the presence of a reflector with solar iron abundance on an infinite plane, 
90% of the plasma has to be blocked from our direct view by an optically thick absorber 
in order to explain the large iron fluorescent line EW measured for V1647 Ori in 2008 by 
Suzaku. Such a geometry is plausible if most of the plasma is hidden behind the stellar core 
and/or disk (see figure 8 and the discussion section in Hamaguchi et al. 2005). We note that 
Aspin et al. (2009) demonstrated a change in the r' band morphology of the reflection nebula 
associated with V1647 Ori (McNeil’s Nebula) between 2004 and 2008, and suggested that 
dust obscuration close to VI 647 Ori may play an important role in defining the observed 
morphology of the nebula. The circumstellar gas structure very close to the stellar core may 
have changed between 2004 and 2008, and this change may have altered both the reflection 
paths of X-ray emission and the openings in visible light between the star and McNeil’s 
nebula. 


We greatly appreciate the Suzaku review committees positive response to our DDT 
request. This work is performed while K.H. was supported by the NASA’s Astrobiology 
Institute (RTOP 344-53-51) to the Goddard Center for Astrobiology. J.K.’s research on X- 
rays from erupting YSOs is supported by NASA/GSFC XMM-Newton Guest Observer grant 
NNX09AC11G to RIT. This research has made use of data obtained from the High Energy 
Astrophysics Science Archive Research Center (HEASARC), provided by NASA’s Goddard 
Space Flight Center. 
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Table 1. Spectral Fits 


Phase 


kT 

z 

N H 

F<i.4 


log Lx 

a 

\ 2 (d.o.f.) 



(keV) 

(solar) 

(10" cm~ 2 ) 

( 10 0 ph cm” 2 s” 

9 (eV) 

(ergs s" 

9 


whole 

4.1 

(3.2— 5.6) 

0.50 (0.31-0.73) 

3.6 (2.9— 4.4) 

2.7 (1.6— 3.8) 

588 

31.1/30.9/ 

'31.3 

1.03 (85) 

early 

2.0 

(1.4— 3.0) 

0.31 (<0.86) 

4.9 (3.4 — 6.6) 

2.4 (0.89-3.8) 

1375 

31.2/30.7/ 

'31.3 

0.65 (36) 

late 

5.3 

(3.7-7. 7) 

0.81 (0.44-1.27) 

3.5 (2.8— 4.4) 

4.6 (2.1 —6.8} 

570 

31.2/31.2/ 

'31.5 

1.25 (47) 

Note. - 

- “A 

bsorption corrected luminosity 

between 0.5— 2 .8 / 

''2.8— 8/0.5— 8 keV, a: 

ss timing d -- 

400 pc. The | 

mrenthe 

ses between 


the second and fifth columns show 90% confidence ranges. 
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Fig. 1. — Left: XIS 0+1+3 color image of the V1647 Ori field (red: 0.5—2 keV, green: 
2—4 keV, blue: 4—10 keV) after correcting the absolute coordinates using a Chandra result 
(Simon et al. 2004). The image was smoothed by the Gaussian function with a =3 image 
pixels (~25") and displayed in a linear scale for the 5 counts pixel -1 range in each color, 
offset at the background level determined from a source free region. The extended blue 
areas at the 3 detector corners are 00 Fe calibration sources. Green circles with numbers and 
a black circle at the center (VI 647 Ori) are X-ray source positions detected in the Chandra 
observation on 2002 November 14 (Simon et al. 2004: Kastner et al. 2004). Radii of these 
circles are proportional to Chandra photon count rates between 0.5—10 keV in logarithmic 
scale, as shown at the bottom right corner. The grey dots show the Chandra ACIS-S FOV. 
Right : monochromatic XIS 0+1+3 unsmoothed image of the same field between 0.5 — 10 keV. 
The cyan and red circles show source and background regions, respectively; the region within 
the magenta circle was excluded from these source and background regions. 





0 2xl0 4 4xl0 4 6xl0 4 8xl0 4 


Time (s) 

Fig. 2. — Top panel: Background subtracted light curve of V1647 Ori between 1-8 keV. 
produced from the XIS 0+1+3 data. Each bin has 2 ksec. Bottom, panel: The hardness 
ratio curve with 4 ksec bins defined by (H - S)/(H + S ), where H and S are the hard 
(4-8 keV) and soft (1-4 keV) band count rates, respectively. The label above the top panel 
shows intervals defined for phase-resolved spectral analysis. The origin of the x-axis is 2008 
October 8. 14 , '33 m 10 4 in UT. 
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Fig. 3.— Left Top: XIS FI (black) and BI (red) spectra of V1647 Ori. The solid lines 
show the best-fit result of the spectra by a model of 1-temperature thin-thermal plasma 
emission (APEC) with a Gaussian function at 6.4 keV for the iron fluorescent line, suffering 
photoelectric absorption by quasi-neutral gas. Left Bottom: residuals of the y 2 values from 
the best-fit model. Right: XIS 0+1+3 spectrum near the iron K line region. The vertical 
axis shows the ratio of the emission against the continuum in the best-fit model (“whole” in 
Table 1). Important emission lines are marked by lines at their rest energies. 



Fig. 4. — X-ray spectra of the early (left) and late (right) phases. 






